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2 : ABSTRACT 

~^ ■ Hydrogen recombination is one of the most important atomic pro- 

cesses in many astrophysical objects such as Type II supernova (SN II) 
atmospheres, the high rcdshift universe during the cosmological re- 
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combination era, and H II regions in the interstellar medium. Accu- 
00 . 
lO ' rate predictions of the ionization fraction can be quite different from 

those given by a simple solution if one takes into account many angu- 

^D | lar momentum sub-states, non-resonant processes, and calculates the 



rates of all atomic processes from the solution of the radiative transfer 
equation instead of using a Planck function under the assumption of 
/\ ' thermal equilibrium. We use the general purpose model atmosphere 

C$ t code PHOENIX ID to compare how the fundamental probabilities such 

as the photo-ionization probability, the escape probability, and the 
collisional de-excitation probability are affected by the presence of 
other metals in the environment, multiple angular momentum sub- 
states, and non-resonant processes. Our comparisons are based on a 
model of SN 1999em, a SNe Type II, 20 days after its explosion. 
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1 INTRODUCTION 

Hydrogen recombination plays a very important role in many astrophysical phenom- 
ena such as Type II supernovae, the interstellar medium, and cosmic recombination. 
Time- dep endent recombination has been stu died for the case of cosmological re- 



combination ( 



2005 



Zeldovich et al. 



Dessart fe Hillier 



1969 



20081 : 



Peebles 



1968 



) and for sup ernovae (jUtrobin &: Chugai 



De. Baron fc Hauschildt 



20091 ). In supernovae this oc- 



curs when the hydrogen recombination time-scale becomes comparable to the age 
of the supernova. Thi s effect is foun d to be dominant in the early epochs of the 
supernova's evolution (jDe et al.l 120091 1 . At later times, however, effects due to multi- 
level atom effects (the importance of having many angul ar momentum s ub-states) 



become more important than time-dependent phenomena. 



De et al 



(|2009|) discussed 



how the effective recombination timescale can be different based on t ime- d ep endent 



De et al. 



rate e quations using different hydrogen atom models. The primary goal of 
( 120091 ) was to determine the epoch in the lifetime of a supernova (during the pho- 
tospheric phase) where time-dependence in the rate equations is most important. 
In doing so we found that at later times model atoms with significantly more en- 
ergy levels (that is additional angular- momentum sub-states) have a strong effect in 
determining the effective recombination time-scale. This issue is also important for 
applications other than supernovae. In fact, considering more complete atomic mod- 
els is important to correctly estimate the electron density and recover subtle features 
in the spectra. Here, we study multi-level atomic systems for hydrogen alone using a 
non-LTE treatment which could alter the hydrogen ionization fraction and therefore 
produce a different temperature structure. 



Cosmological recombination codes such as RECFAST (ISeager et al. 



1999 



20001 : 
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Wong et al.ll2008l ) and RICO (IFendt et al.ll2009l ) that deal with cosmological recombi- 
nation and solve for the free electron fraction as a function of redshift do an excellent 
job. Nevertheless there a re assumptions about Saha equilibrium between higher an- 



gular momentum states. 



Grin &: Hiratal (120 lOl ) have presented their recombination 



code RECSPARSE which is a highly detailed calculation but still uses the Sobolev ap- 
proximation for line calculations. Recently the cosmological recombination problem 
has been revisted by several authors, paying special attention to possibly neglected 



atomic p hysics effects. The processes studied have inc 



and He I ( 



Hirata 



Dubrovich fc Grachev 



2008; 



2005 



Karshenboim fc Ivanov 



Wong fe Scott 



2008 



absorption of the He I A584 line photons (ISwitzer fc Hirat 



uded two-photon decays in H 



2007 



Chluba fc Sunvaev 



Hirata fc Switzer 



2008 



Kholupenko et a 



sorption ( 



Hirata 



man scattering ( jHirata 



2008 



al hoOSai bl: 



2008 



2008 ). H I continuum 



Rubino-Martfn et al. 



20071). stimulated two-photon decays and two-photo n ab 



Chluba fc Sunvaev 



2008 



gular momentum su 



( Wong fc Scott 



jstates ( 



2007 



2006 



Hirata fc Switzer 



Grin fc Hirata 



Switzer fc Hirata 2008a. 



20 



Karshenboim fc Ivanov 



2008 



20081 ). Ra- 



a more complete set of an- 



01), and forbidden transitio n s in He I 



Rubino-Martin et al. 



2008|). 



Hirata fc Forbes 



( 120091 ) studied the effects of partial frequency redistribution, line blanketing, and 
time-dependence in Ly-a, but even then made use of the Sobolev approximation. 
Here we resolve all line profiles in the co-moving frame and accurately treat line 
blanketing. We assume complete redistribution for the line profiles. 

We focus primarily on estimating the photo-ionization rate, the escape probabil- 
ity, and the collisional de-excitation rate. We show how these rates depend on the 
hydrogen model atom used and we also study the effects of the presence of metals 
on the results. These different rates ultimately lead to variations in the hydrogen 
ionization fraction. We also examine when the 2j process plays the most significant 
role in controlling hydrogen recombination. 

In § [21 we outline the theoretical framework relating the photo-ionization rate, 
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collisional de-excitation rate, and escape probability to the effective hydrogen recom- 
bination time. In § [3j we describe our approach and different test systems that we 
study In § H] we present our results. In § [5J we state our conclusions and describe a 
possible framework for future work. 

2 THEORETICAL FRAMEWORK FOR HYDROGEN 
RECOMBINATION 

The ionization fraction for hydrogen is defined as 

/ H = ^S_. (1) 

nm + n H 

For any system the ionization fraction depends on the net recombination rate in 
the system. For a system in LTE, the ionization fraction can be calculated using 
the Saha equation with knowledge of the electron density and temperature. On the 
other hand, for systems not in LTE, it is important to solve the radiative transfer 
equation and the rate equations simultaneously. The recombination rate of a system 
is dependent on the number of angular momentum sub-states of the system and also 
on the presence of other elements in the environment that can contribute to the free 
electron density. The recombination process in a system is determined by the net 
rate of photo-ionization and its inverse process (recombination) which are modified 
due to line transitions. The effectiveness (towards recombination) of a resonant tran- 
sition depends on the escape probability of the line. Also there are non-radiative 
transitions that take place such as the collisional de-excitation of electrons to lower 
energy levels. In addition to these there is the downward 27 process that connects 
states of the same parity while the corresponding upward transition rate is very low. 
Therefore it is important to study all these quantities to determine the correct na- 
ture of recombination. We investigate the 27 process, the photo-ionization rate, the 
escape probability and the collisional de-excitation rates of systems with different 
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hydrogen atom models and metallicities of the environment. We denote Pn as the 
photo-ionization rate from a state with index n (which is not its principal quantum 
number). Similarly, we define the escape probability as /3„i which is the escape prob- 
ability for the resonant line photon between the level Is and the level characterized 
by index n. The collisional de-excitation rate is defined as the rate of de-excitation 
of the electron from level n into the Is state. The analytical evaluation of each of 
these quantities is possible if the system is in LTE or if the system has a model atom 
which has just a very few angular momentum sub-states. But for a system with many 
angular momentum sub states (we consider 921 bound levels for hydrogen) and in 
the presence of metals, it is necessary to evaluate these quantities numerically from 
a consistent solution of the radiative transfer equation and the rate equations. 

The results presented in this paper are obtained by simultaneously solving the 
coupled radiative transfer equation and rate equations. The rate equations use multi- 
level atoms including a large set of line and continuum transitions. The upward 
radiative rate (i —¥ j) is given as niRij where i is a bound state and j can be a 



boun d or continuum state, n, is the population for level i. Rij is given as (iMihalas 



1978|): 



Ra = Air 



a.ij(v)(hv) 1 J v du 



(2) 



Similarly the downward rate (j —¥ i) is calculated as nj(—)*Rji where Rji is: 

,2/uA 



R 



j i 



4-7T 



"o 



aij(u)(hu) x 



+ J V 



his 

e kT dv . 



(3) 



In the above equations J„ is the angle averaged radiation intensity at frequency v. T is 
the temperature, k is Boltzmann's constant, a^ is the photo-ionization cross section 
for the i — )- j transition, and uq is the threshold frequenc y. The level pop ulations 
denoted with (*) are the equilibrium values as defined by iMihalasI (119781 ) . In our 
calculations the Sobolev approximation has not been used. In order to compare our 
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results with others we plot the values of the Sobolev escape probability. All the level 
populations were calculated without assuming the S obolev approx imation. 



Below we give the equations that define P n , f3 n i (IMihalas 



19781 ). Firstly, 



ft 

n K 
where k stands for the continuum level. 

The escape probability f3 n \ is defined using the Sobolev escape probability as 

0m = , (5) 

Vnl 

Vm = fniKiitlni -n n ) , (6) 

m e c \ g n ) 

where f n \ is the oscillator strength of the line between excited level n and the ground 
state, ni and n n are the number densities of the ground state and the excited state 
referred by index n, respectively. A n i is the wavelength of the line, t i s the time sinc e 



explosion in seconds. g n is the degeneracy of level n and — ~ -\ ( JMihalaslll978l ). 



In § [3] we discuss our approach to quantify the photo-ionization rate, the escape 
probability, and the collisional de-excitation rate. 

3 VARIATION OF HYDROGEN MODEL ATOMS AND 
COMPOSITION 

3.1 Description of the test systems 

To set up the physical structure of each system we use the density profile and lumi- 
nosity from a full NLTE calculation of a model atmosphere in homologous expansion, 
with a power-law density profile p oc v~ 7 , 20 days after explosion, which is a rea- 
sonable fit to observed spectra of SN 1999em at that epoch. This underlying density 
profile and the total bolometric luminosity in the observer's frame are chosen to be 
representative of the conditions in SNe IIP near maximum light. Give this structure 
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we use our general purpose model atmosphere code PHOENIX ID to solve for the 
new temperature structure and level populations under the new conditions (such as 
the hydrogen model atom and the presence of other metals). The initial model of 
SN 1999em at 20 days after explosion was generated using a 31 level hydrogen model 
atom and solar metals were included. We call this starting model our base model. 
Our base model was generated treating hydrogen and ot her elements in NLTE. The 



details of our radiative transfer code are described in 



Hauschildt fc Baron 



( 11999 



20041 and references therein). The transfer equation and the rate equations were all 
time independent for the base model. The physical properties related to our base 
model are given in Table [U 

Now we define the different systems that we experimented upon to study the 
behavior of the physical quantities affecting the recombination rate. For convenience, 
we name the different systems in the following way. We refer to Model A as the 
system which has only hydrogen and the hydrogen model atom in this case has just 
4 bound states (Is, 2s, 2pi,2ps) and 4 lines that couple those states. Model C is 
the same as Model A, except it has a solar composition of metals present in the 
system in addition to hydrogen and all the other metals are treated in LTE. Model 
B has only hydrogen in the system but the hydrogen model atom has 921 bound 
levels and 995 lines that couple the bound states. Model D is the same as Model B 
except this has metals present and the metals are treated in LTE as in Model C. The 
LTE treatment of metals was chosen to reduce computing time and in a hydrogen 
rich environment the main role that metals play in continuum processes is to be a 
source of additional electrons. We did not focus on how their atomic structure affects 
hydrogen recombination. 

For Models A and C, there are primarily three major processes considered that 
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couple the bound states: 

His + 7 -»• e~ + H+ ; 
H ls + j <-► H 2p ; 

H 2s ->• #is + 2 7 . 

In addition to this the Is and 2p states are coupled b y spontaneous transitions. In 



models B and D the atomic line data were taken from 



Kuruczl (Il995[ ). 



4 RESULTS 

We study the rates and the resulting hydrogen ionization fraction, f H , among the 
four different systems that have been described in § [31 These different systems can 
have different ionization fractions primarily because the solution of the transfer equa- 
tion and the resulting temperature structure will be different than in the base model. 
The models that have metals (Models C and D) are expected to have temperature 
structures similar to the base model. The 4-level hydrogen atom model without met- 
als is likely to have a different temperature structure than the base model. These 
differences could be subtle but important enough to change to the hydrogen ion- 
ization fraction. In our results we have also studied whether the 2 7 process has a 
significant effect on the recombination time scale or the resulting ionization level or 
the consequent temperature structure of the system. 

We define T st d as the optical depth corresponding to the continuum opacity at 
500 nm. We classify our results depending on whether the temperature structure was 
held fixed at the base model value or was allowed to reach radiative equilibrium under 
the particular assumptions of each case, in order to isolate the effect of temperature 
change on the level populations and the hydrogen ionization fraction. 
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Table 1. Physical characteristics from our base model. 

t p (g cm~ 3 ) R (cm) v (km s _1 ) T (K) n e (cm -3 ) p 

9.897E-05 1.823E-15 2.454E+15 14200.00 3157.61 99464.3 1.34 

3.626 3.084E-13 1.242E+15 7187.50 7345.20 1.079E+11 0.75 

804.480 6.533E-12 4.104E+14 2375.00 26155.8 2.97E+12 0.67 

The columns give various quantities as a function of optical depth, p is the mass density, R is the radius, v is the 
velocity, T is the matter temperature, n e is the electron density, and p is the mean molecular weight. 

4.1 Metal-rich Models 

We first discuss our results for the models where the composition included metals. 
We treat helium and metals in LTE. Hydrogen is always treated in NLTE. Recall that 
these systems are termed Models C and D for the 4- level and 921-level hydrogen atom 
models, respectively. Figure [T] shows the hydrogen ionization fraction for Models C 
and D. There is a significant change in the hydrogen ionization level, /#, between 
Models C and D. The quantity fn decreases in the multi- level atom case in the lower 
optical depth regime. For r st d > 0.1 the ionization levels among Models C and D 
are not very different. The quantitative difference is also tabulated in Table |2] which 
shows the physical parameters for Models C and D when the 2j process was or was 
not included. The reduction in /# due to additional angular momentum sub-states 
was about a factor of 3 at an optical depth of about r st d ~ 10~ 4 (when the 27 process 
was included in both the models). The difference in /# decreases as the optical depth 
increases. In metal-rich systems, (Models C and D), the exclusion of the 27 process 
did not seem to have any significant effect for almost all optical depths of interest. 

Figure |2] displays the net photo-ionization rate from any bound state of the hydro- 
gen atom, P n , for the levels 2pi/ 2 and 2p 3 / 2 for Model C The profile of P n for Model 
C is not monotonic, but there is an overall trend to increase with optical depth, r s td- 
This increase in P n is expected due to the fact that photo-ionization dominates over 
recombination at higher optical depths due to higher temperatures. Figure [3] shows 
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the net photo-ionization rate as a function of wave number of each energy level for 
Model D. Different panels show different T st d regimes. The P n values increase with 
increasing r st d for a given energy level. Also the net photo-ionization rate does not 
change significantly with the change in the energy level of the bound state. There is 
a drop in the P n profile at energy levels very close to the continuum. This may be 
due to the fact that very high energy bound states are not really distinguished from 
the continuum, therefore the net photo-ionization rate falls off. 

For the calculation of the escape probability (defined as the probability for the 
escape of the resonant line connecting the ground state and higher energy bound 
state), are shown in Figures HH5J Figure H] shows the escape probability for Model 
C, for levels 2pi/ 2 and 2p 3 / 2 . The escape probability for Model C decreases with 
increasing optical depth and also decreases with decreasing energy of the bound state. 
For Model D (Figure [5]) the escape probability increases with increasing energy of 
the bound state. Also the escape probability decreases very slowly with the increase 
in optical depth for Model D, similar to that in Model C. The apparent difference 
between Models C and D is that the escape probability decreased for a higher energy 
level in Model C as opposed to case D where the escape probability increased with 
the increased energy of the bound state. For Model C, the levels in question are both 
2p states and in Model D for those two 2p states we also find a slight reduction in 
the escape probability with increasing energy of the bound state. 

The reduction in the escape probability as a function of optical depth can be 
explained from Eq. EJ At higher optical depths (for a given excited state), when n n , 
the population of the excited state increases, this results in a decrease in the quantity 
rj nl . Following the decrease in r) n i, the escape probability f3 n \ decreases. On the other 
hand at a given optical depth with increasing energy of the bound state, the escape 
probability increases. At a given r st d-, for an increase in the bound state energy, n n 
decreases and A n i decreases, hence the overall effect almost always is to decrease (5 n \. 
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Figure |6] shows the collisional de-excitation rate (which is the collisional de- 
excitation coefficient q e times the free electron density n e ) for Model C. In this Model 
there is an increase in the collisional de-excitation rate with increasing optical depth. 
For Model D, the collisional de-excitation rate follows a similar pattern. The increase 
in the free electron density with increasing r s td explains the increased collisional de- 
excitation rate with r st d- To summarize, the characteristics of the metal-rich Models 
C and D, 1) the ionization fraction decreases for a multi-level hydrogen atom com- 
pared to a 4-level hydrogen atom. 2) The net photo-ionization rate increases with 
optical depth for a given bound state. The photo-ionization rate also does not change 
significantly at a given r s td as the energy of the bound state increases. 3) The escape 
probability decreases with increasing T st d, although this decrease is very small. The 
escape probability increases with increasing energy of the bound state except for the 
case of the 2p states. 4) The collisional rate increases with increasing optical depth. 
Also for the multi-level case the rate is close to the equilibrium value. 5) The 27 
process did not seem to have any significant effects in Models C and D. 



4.2 Metal-deficient Models 

In this section we focus on Models A and B. Recall that Models A and B are models 
where the atmosphere consists of pure hydrogen. Hydrogen was treated in NLTE in 
these models. In Model A, we use a 4-level hydrogen atom and in Model B, we use a 
hydrogen atom model with 921 energy levels. We also study how the 27 process affects 
these systems when the environment is metal-free. Figure [7] displays the ionization 
fraction for Models A and B. Interestingly fjj does not change significantly between 
Models A and B. This is quite different when compared to Models C and D. Table [2] 
shows the changes due to the 27 process. For Model A, at low optical depths (r st d ~ 
10 -4 ) turning on the 27 process produces about a 20% change in the free electron 
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density, but for r st d > 10~ 3 no significant change is seen. For Model B, there is a 
significant reduction in the free electron density due to the inclusion of the 27 process. 
The electron density is 5 times higher with the inclusion of the 27 process in case B, 
at low optical depths, r st d ~ 10~ 4 . For higher optical depths, there is no significant 
effect due to the 27 process. The decrease in the free electron density due to the 27 
process is expected since the upward process is suppressed by a factor of 1/137. The 
change is especially noticeable in Model B. It is at a fairly low optical depth where 
the change in the rates observed but this is also the optical depth where the optical 



to what is described in 



De et al. 



depth of the Balmer li ne is high ( De et al.l 120091 ) . This wavelength regime is similar 



f l2009h 



The net photo-ionization rate (P n ) does not change significantly between different 
bound states in both Models A and B (see Figures [8] and [9]) . Although in the multi- 
level case (case B), P n increases by almost a factor of 10 from the lowest energy 
bound state to the higher energy bound states (except for the states very close to 
the continuum) at low optical depths (r s td < 0.01). Figures [10] and [TT] show the 
escape probability for Models A and B respectively. The trend is very similar to 
that of Models C and D. The collisional de-excitation rate increases with increasing 
optical depths in both the Models C and D 

To summarize: our findings on Models A and B (pure hydrogen models) 1) For 
most optical depth regimes, the basic trend in the rates is similar to Models C and 
D. 2) The 27 process seems to have a significant effect in both Model A and B (at 
lower r st d)- The effect is much larger for Model B. This effect was not seen in Models 
C and D. 

We observe the following, for almost all optical depths, 

P n {A) > P n (C); 
f3 nl (A) < p nl (C); 
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n e q nl (A) < n e q nl (C). 

For Models B and D, 

P n {B) < P n {D) (r<0.01); 
MB) « /3 nl (£>); 
n e q nl (B) < n e q nl (D) (r < 0.01). 

The difference in the profiles of P n , /3 n i and n e g n i between Model A (or B) and Model 
C (or D) can be attributed to their different temperature profiles. 

4.3 Metal-deficient case without temperature corrections 

In our previous Models A-D, we allowed the system to reach radiative equilibrium 
for each assumption. It is useful to hold the temperature structure fixed and just 
examine the changes that are due to the variation in the compositions and model 
atoms. In the fixed temperature structure case we do not find a large difference 
as compared to the radiative equilibrium case. The differences, compared to the 
radiative equilibrium case occurs in the optical depth region 10~ 3 < r < 10~ 2 . There 
is a slight suppression in the ionization fraction for Model B compared to Model A. 
The other rates are nearly identical. 

5 DISCUSSION 

The primary motivation for this paper was to investigate how the multi-level struc- 
ture of the hydrogen model atom affects the hydrogen ionization fraction as well as 
the photo-ionization, escape, and collisional de-excitation rates. We also investigated 
if the the metals in the environment or if the 27 process plays a significant role in 
the recombination process of the system. We find that a multi-level hydrogen atom 
structure, or in other words a more complete set of energy levels to be important 
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Table 2. Physical characteristics from four different models, under radiative equilibrium 



Model Name 



Type 



T (K) JV e (cm" 3 ) 



4 level 



0.00016 

0.02343 

3.61 

1104.90 



2820.75 

4011.17 

7564.24 

27565.10 



105236. 
9.999E+07 
2.351E+11 
3.903E+12 



1.01 
1.01 
0.54 
0.50 



4 level, no 27 process 



0.00016 

0.02345 

3.61 

1106.10 



2820.61 
4011.18 
7564.29 
26840.2 



121450. 

9.999E+07 

2.351E+11 

3.9035E+12 



1.01 
1.01 
0.54 
0.50 



921 level 



0.00016 

0.0272 

1.794 

1181. 



2801.08 

4547.80 

6779.06 

20462.42 



3.946E+5 

2.981E+8 

1.087E+11 

3.893E+12 



1.01 
1.01 
0.57 
0.50 



921 level, no 27 process 



0.000157 

0.02705 

1.794 

1181 



2826.28 

4547.69 

6778.83 

20462.42 



6.709E+4 

2.983E+08 

1.0864E+11 

3.801E+12 



1.01 
1.01 
0.57 
0.50 









0.00023 


3564.79 


271828. 


1.26 


c 


4 level, 


metal-rich environment 


4.91280 


7584.85 


1.260E+11 


0.68 








867.030 


26365.40 


3.162E+12 


0.625 








0.00023 


3564.78 


271808. 


1.26 


c 


4 level, 


metal rich environment 


4.91270 


7584.84 


1.260E+11 


0.68 






no 27 process 


867.900 


26243.10 


3.160E+12 


0.625 








0.00011 


3099.54 


84640.4 


1.26 


D 


921 level, 


metal rich environment 


0.98338 


6337.20 


4.254E+10 


0.98 








863.240 


26109.80 


3.156E+12 


0.63 








0.00011 


3099.46 


84638.0 


1.26 


D" 


921 level, 


metal rich environment 


0.98373 


6337.44 


4.256E+10 


0.98 






no 27 process 


863.080 


26113.50 


3.156E+12 


0.63 



in determining the ionization profile of the system. Recently, iGrin fc Hiratal (120101 ) 
studied the effects of including a more complete set of angular momentum substates 
in the study of cosmological recombination and also found that including more sub- 
states to be important. The need to have a more complete set of bound states in the 
hydrogen atom is found to be more important in a metal- rich environment. This is 
because of the fact that recombination is more effective in a multi-level framework. 
When there is a source of additional free electrons (such as metals), the suppression 
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in Jh is larger than in pure hydrogen models. This is merely due to the fact that 
the larger free electron density drives recombination, thus it is very important to use 
multi-level hydrogen atoms especially in realistic environments with solar composi- 
tions. 

The enhancement of recombination is indicated by the increasing escape proba- 
bility for the higher bound states (for both Models B and D). This is also reinforced 
by the almost constant photo-ionization rate over different energy bound states. 

The importance of the 27 process is seen in the 4-level pure hydrogen case (Model 
A) as only a small reduction in the free electron density. In Model B, the reduction in 
the free electron density due to the inclusion of the 27 process is much larger. There 
was not any significant effect due to the 27 process in the solar composition Models. 
The larger change in Model B (at lower r st d) can be explained from Figure [121 At low 
optical depths the relative population of the ground state to the first excited state 
is displayed. For the metal-free Models (A and B) the ratio — is not very high and 
about 4.0 for Model B. Thus in Model B, the first excited state population is not 
significantly lower than the ground state population. The relative first excited state 
population is much lower in Models C and D (where the ratio is around 10 4 — 10 5 
at similar T s td)- The high level population of the first excited state (for Model B) 
increases the importance of the 27 process, due to the higher number of available 
electrons that can recombine into the ground state and the upward transition prob- 
ability for this non-resonant process being low, thus the recombination obtained is 
more effective. We find that it is essential to incorporate both multi-level atoms 
and the 27 process. Our primary purpose was to do a simple calculation to empha- 
size the important factors in determining the recombination of hydrogen at typical 
free electron densities found in Type II supernova atmospheres. Our results show 
that it is important to incorporate 27 transitions and multiple angular momentum 
sub-states at low optical depths where the free electron density is small (typically 
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~ 10 5 cm" 3 ) and is also ty pical of the fr ee electron density during the epoch of 



cosmological recombination ([Peebles 



19681 ). In H II regions the electron density is 



much lower making the hydrogen recombination time scale rise. Therefore accurate 
treatment of hydrogen recombination is important in all of these scenarios. In the 
future we plan to investigate hydrogen recombination in the context of cosmological 
recombination epoch. 
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- 1 0.10 - 



Figure 1. Comparison of the hydrogen ionization fraction found using the 4 level (C) and 921-level (D) model atoms 
in a metal-rich environment. The upper panel shows the lower optical depth regime while the lower panel shows the 
higher optical depth regime. We define T stl j as the optical depth corresponding to the continuum opacity at 500 nm. 
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Figure 2. Photo-ionization rates of the 2p states versus optical depth for Model C. 



© 2010 RAS, MNRAS 000, 000-000 



20 De, Baron, & Hauschildt 




Figure 3. Photo-ionization rates versus energy level for Model D at different optical depths. Each panel shows a 
particular optical depth. 
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Figure 4. Escape probability versus optical depth for the 2p states for Model C. 
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Figure 5. Escape probability versus energy levels for Model D at different optical depths. Each panel shows a 
particular optical depth. 
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Figure 6. Collisional dc-excitation coefficients of the 2p states versus optical depth for Model C. 
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Figure 7. Comparison of hydrogen ionization fraction obtained using the 4 level (A) and the 921-level model (B) 
hydrogen atom in pure hydrogen under radiative equilibrium. The upper panel shows the lower optical depth regime 
while the lower panel shows the higher optical depth regime. 
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Figure 8. Photo-ionization rate versus optical depth for 2p states of Model A under radiative equilibrium. The last 
point at r 3 td = 10 is a numerical glitch caused by poor spatial resolution. The rate should continue to drop with 
depth. 
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Figure 9. Photo-ionization rates versus energy level for Model B in radiative equilibrium. Each panel shows a 
particular optical depth. 
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Figure 10. Escape probability versus optical depth for the 2p states of Model A in radiative equilibrium. 
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Figure 11. Escape probability versus energy level for Model B in radiative equilibrium. Each panel refers to a specific 
optical depth. 
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Figure 12. Ratio of the ground state to the first excited state population densities as a function of optical depth. 
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